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Abstract: We perform a parton level study of angular correlations, in higher dimensional 
black hole scenarios at the LHC. We start by discussing some features of the angular 
Q-i spectrum in the high energy limit using the geometrical optics approximation. This allows 

c~| , us to recover the high energy limit of the Hawking fluxes. Then we use the full Hawking 

angular fluxes on the brane, for a singly rotating Myers-Perry black hole, to motivate 
the construction of various angular observables as to maximise the angular asymmetries 
due to rotation. This is finally adapted to a parton level simulation using the CHARYBDIS2 
generator. We explore two types of variables based on: axis reconstruction and two-particle 
angular correlators. We find energy cuts which have the potential to help identify the effects 
of rotation in semi-classical rotating black hole events at the LHC during the 14 TeV run, 
especially for spin-1 particles. 
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1 Introduction 

Theories with extra dimensions have been proposed to address the hierarchy problem of 
the standard model of particle physics (SM) [1-6]. In such theories, the scale of gravity 
is lowered to the TeV energy scale and strong gravity phenomenon becomes possible, in 
particular black hole production as suggested by several authors [7—11]. Theoretically, this 
is further supported by evidence for black hole formation in higher dimensions from several 
analytic studies [11-16] and numerical relativity studies [17-21]. 

From the point of view of detection, such states are predicted to decay instantaneously 
through Hawking evaporation [22-32] if produced at the Large Hadron Collider (LHC). 
Though there are still theoretical uncertainties such as the amount of gravitational radi- 
ation emitted during the collision [15, 16] and the corrections of the Hawking spectrum 
during evaporation due to gauge charges on the brane [33, 34] or the mass of the parti- 
cles [33-35], current event generators already take into account the full angular spectrum 
with black hole rotation. Two such event generators, CHARYBDIS2 [36] and BlackMax [37] 
are currently being used at the LHC to look for signatures of black hole production and 
evaporation. Bounds on these models, particularly on the fundamental planck scale Mu 
have emerged over the years, up to the present LHC data. The main sources of bounds 
come from: i) deviations from Newtonian gravity in torsion balance experiments [38, 39], 
ii) collider searches for Kaluza-Klein graviton production (monojet or photon with missing 
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energy) [40-42] and KK graviton mediated dilepton, diphoton or dijet production [43-48] 
iii) astrophysical or cosmological KK graviton production in supernovae and in the early 
universe [49, 50]. The most restrictive laboratory bounds from current LHC searches in- 
dicate Mo > 2.5 — 3.5 TeV for D > 6 whereas observational bounds from astrophysics 
and cosmology only allow for such small for D > 7. These bounds should be taken 
as indicative because they are very model dependent, especially the astrophysical and cos- 
mological bounds. Regarding direct searches for black hole events at the LHC, the first 
attempt to establish some bounds for these models was released by the CMS [51] and AT- 
LAS [52, 53] collaborations. However, their analysis for the 7 TeV data, depends on regions 
of parameter space where black holes with masses close to the unknown Planckian regime 
would be produced [54]. If cuts are imposed, such that the semi-classical approximation is 
valid, the cross-sections become negligible at 7 TeV, so this channel will only be properly 
tested in higher energy runs (originally planned to be of 14 TeV) which are likely to start 
by the end of 2014. 

In this article, we focus on the semi-classical limit for the 14 TeV LHC runs and explore 
the possibility of observing angular asymmetries in the large extra dimensions scenario. 
We use the default model implemented in the CHARYBDIS2 event generator, which takes 
into account all of the best known theoretical calculations with black hole rotation. The 
purpose of this study is to focus on how to extract the physics from the events at parton 
level, leaving aside a detailed hadron level analysis and detector simulation. 

Our study starts with an analysis of the high energy limit of the Hawking spectrum 
through a geometrical optics calculation. It is well known that in the high energy limit 
wave propagation can be treated geometrically. Therefore, in this limit, the study of the 
classical trajectories of test particles outside the black hole provides useful information. 
The aim is to obtain the shapes of the absorptive discs as seen by an observer away from 
the black hole 1 This calculation provides: i) some physical intuition which helps visualising 
the background and ii) the correct high energy limit of the Hawking spectrum, which could 
be useful to implement high energy tails in current event generators. 

Next, we provide a summary of the full spectrum of Hawking radiation on the brane, 
for a singly rotating Myers-Perry Black hole which is well known in the literature [22, 24- 
27, 30-32, 56]. This contains crucial information which guides our construction of various 
angular variables. The full parton level modelling of the events is done in the second 
part using the CHARYBDIS2 event generator which includes, in full, the effect of rotation 
for brane degrees of freedom. We study two types of observables. The first type is an 
attempt to reconstruct the angular momentum axis of the collision. Guided by the angular 
distributions with the true angular momentum axis, we study a reconstruction based on the 
hardest particle and the sphericity tensor of the event. Though there is some correlation 
with the true axis we find it fails for a considerable fraction of events. The second type 
of variable is much more promising. We use angular correlators between pairs of particles 
with definite helicity. This explores the strong helicity dependence of the Hawking angular 

X A first study in higher dimensions was done for some special cases in [55]. We generalise their arguments 
to an arbitrary trajectory. 
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fluxes and its dependence with energy, especially at small energies (where particles with 
opposite helicities tend to be emitted in opposite hemispheres). We find that with suitable 
energy cuts, we can obtain a strong asymmetry. We suggest that the effect should survive 
when a full simulation is done if we explore the helicity dependence of W and Z decays, or 
use low energy jets. 

The structure of the paper is the following: In section 2 we present the background 
spacetime and present the geometrical optics limit calculation. We discuss how such a 
calculation matches well with the full Hawking spectrum and present the full Hawking 
spectra implemented in CHARYBDIS2. In section 3 we construct several angular variables 
presenting an analysis of some Monte Carlo samples which show how to extract the effects 
of rotation at parton level. We summarise our conclusions in section 4. Some details of 
the calculations are left to the appendices. 

2 The Background 

In this article we are interested in the effects of rotation on the angular emission spectrum. 
Though in general a brane world black hole may contain additional gauge charges (for a 
study see [33, 34]) we consider the simplest case of a singly rotating Myers- Perry black hole 
formed on the brane with line element (D = 4 + n) 

r 2 + a 2 -A\ 2 (r 2 + a 2 - A) 2a sin 2 9 £ l2 



ds = ( 1 - ) dr + ^ dtd(j) - —dr 



where 



n / o n at (v 2 + a} — A) Sffi 2 9 \ 

- Y,d9 2 - \r 2 + a 2 + — i '- sin 2 9d(\> 2 - r 2 cos 2 9dQ 2 n , (2.1) 



A = r 2 + a 2 - , S = r 2 + a 2 cos 2 9 . (2.2) 



The parameters {/i, a} are related to the physical mass M and angular momentum J given 
by 

M (n + 2) n(n + l) . i 

- ^i5 2+n (27r)-— M™+V , (2.3) 



J = 5 2+n (2vr)^^+^ Ml +2 aji = —Ma , (2.4) 

n + z 

where Md is the higher dimensional Planck mass and S'2+n is the area of the n + 2-sphere. 
The horizon radius and oblateness {rjj, o*} are defined through the largest positive root of 
A(r#) = (rjj is directly related to the surface curvature of the horizon) and a* = a/rn 
(which is the oblateness of the spheroidal horizon) respectively. 

The theory of Hawking radiation predicts that black holes emit a continuous flux of 
particles. In horizon radius units (rjj = 1), the differential fluxes of particle number, energy 
and angular momentum are [57] 

d{N,E,J} s 1 ^ A {l,w,m} m a , m ,2 , n K \ 

r- ; — 77^- = — > > ; — ; ; ; r^rrTIi '{OJ, a) Lb (aco, cos 0) , (2.5) 
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where u> = u — m^ln, k = {j, m} are the angular momentum quantum numbers, h 
the helicity of the particle, 

(n+ 1) + (n- l)a 2 



47r(l + a 2 )r# 



-s is 



(2.6) 



and £Ih = a/(l + a 2 ) is the angular velocity of the horizon. Y^ +n ^ are the so called 
transmission factors defined as the fraction of an incident wave from infinity which is 
absorbed by the black hole. The boundary conditions are such that close to the horizon 
the wave is purely ingoing for co-rotating physical observers [58]. The term in denominator, 
is usually called the Planckian (or thermal) factor. In the next sub-sections we analyse 
some of the properties of the Hawking radiation that are expected in this background. 

2.1 The high energy limit 

We start with a study of the geometrical optics (or test particle) limit. We focus on 
geodesies for particles on the brane, where the brane metric is obtained through a projection 
dfi 2 = 0. A classical particle stuck on the brane follows a geodesic curve x M (A), with 
parameter A, determined by varying the action 



S 



1 dx a dx a 
1 2 ~dX ~dX 



(2.7) 



In this formulation, the conserved quantities are identified by looking at the symmetries 
of the Lagrangian, or equivalently, the Killing vectors of the metric. The brane metric, 
given by the projection of (2.1), has the same form as the Kerr metric. This type of 
metric produces two obvious conserved quantities associated with its time and azimuthal 
Killing vectors and a third one related to a Killing tensor. For example in [55] the three 
conserved quantities were combined with the Hamiltonian (which is also conserved since 
the Lagrangian (2.7) does not depend on A) to obtain a radial equation of motion for a 
particle with mass [i. We can apply the same reasoning to our case, where we will consider 
more general trajectories. For simplicity we switch to horizon radius units where ru = 1 
so that fi — > 1 + a 2 and we have the following mapping of parameters 
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The equations for the geodesies are then 
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Incident particle 




Figure 1. Diagram defining the £ angle and the impact parameter b on the plane transverse to the 
direction of incidence. 



where uj,£ z ,Q are the constants of motion associated with time translations, azimuthal 
translations and the Killing tensor respectively. Next we redefine A — > Acj, express the 
remaining constants of integration in terms of the impact parameter b (in units of r#), the 
polar angle of incidence i9, the angular momentum magnitude £, the angular momentum 
orientation at infinity and the reduced mass v = fi/u : 
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The impact parameter b corresponds to the distance of closest approach to the origin, if the 
spacetime was flat. £ is the polar angle of the impact parameter on the plane perpendicular 
to the direction of incidence (see figure 1) or equivalently the angle between the angular 
momentum of the black hole and the angular momentum of the incident particle. 

The most general trajectory is parametrised by the set {6, £, a, u}. Our goal is to 
determine whether a particle is absorbed or not for a given set. This can be achieved by 
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looking at the radial equation of motion in (2.10), and noting that 1Z must be non-negative 
over the trajectory. Let's start by defining the functions A,B,C through 

n = — [-Ab 2 + 2Bb + C] . (2.12) 

To investigate whether a particle is absorbed or not, we need to know if 1Z takes negative 
values at any point of the trajectory. If so, there is a region which is inaccessible. A way to 
guarantee absorption is to choose the value of b in an interval such that 1Z is positive over 
the whole space. This requires determining its zeros with respect to b at all points r and 
finding the intervals of b such that such zeros do not exist over the whole radial domain. 
In principle this may produce complicated discs or ring- like regions on the {b, (,} polar 
plane at infinity. In what follows we focus on v = 0, though the treatment can be easily 
generalised. This is because v = /i/w, so since the geometrical description is good for u 
large, this parameter should be small. The detailed analysis in appendix A leads to the 
following result. The range of impact parameters for absorption can be summarised as 

< b < b c = min j = 1 (2.13) 

{r>r low } I VWTAC -B) 

with 

ri ow = < 




and ta defined such that A\ rA = 0. 

Our goal is to obtain the absorptive discs as seen from infinity in various directions 
of observation. The most efficient way is to numerically determine the critical impact 
parameter for each set of parameters through the minimisation of (2.13). Another useful 
approach is to solve the problem analytically in a particular case and expand around it 
perturbatively. This method is useful as a check of the numerical analysis and shows 
most of the features of the result. Equation (2.13) is easy to minimise when a = (the 
Schwarzschild case), leading to 



ro = 



n + 3\4i ^=(IL^^\/!E|, (2.15) 



ro, bo are the minimiser and minimum respectively. For small rotation parameters, it should 
be possible to find a good approximate solution of (2.13), by expanding perturbatively 
around (2.15), i.e. 



hoo oo , 



b c EE V (2.16) 



p! ^— ' m! 

p=0 r m=0 

with r c , b c the minimiser and minimum for a / 0, respectively. These expansions can be 
inserted in (2.13) as to find the coefficients order by order in perturbation theory. The 
recurrence relations are provided in appendix A.l. 

Figures 2 and 3 show a perfect agreement between the pertubative method and the 
numerical method, for small a and they are good up to a ~ 1. For larger a, the perturbative 
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n = 2, a = 0.1, 8 = § n = 4, o = 0.5, 9 = J n = 10, a = 0.5, 9 = § 

2 | , , , , , , , 1 2 | , , , , , , , 1 2 | , , , , , ,— 




I i i i i i i i I _2 I i i i i i i i I _2 I i i i i i i i I 

-2 -1.5 -1 -0.5 0.5 1 1.5 2 -2 -1.5 -1 -0.5 0.5 1 1.5 2 -2 -1.5 -1 -0.5 0.5 1 1.5 2 



Figure 2. The perturbativc result (dashed blue) agrees with the numerical minimisation (solid 
red), for a < 1 (some combinations of (n,a,"$) indicated). The polar system (6, C) used throughout 
to draw the absorptive discs (defined on the plane transverse to the radial direction of observation 
at infinity) is indicated on the left plot. 



n = 4, a = 1, = J, 0(6) n = 4, a = 1, 9 = f , 0(2) n = 4, a = 2, 9 = J, 0(2) 




Figure 3. The pcrturbative result starts to fail compared to the numerical minimisation, for a > 1 
(same color scheme as figure 2). n = 4, a = 1 (left) starts to disagree if we include 0(6) corrections. 
However, truncating at 0(2) (centre) the result is still good. Increasing a = 2 (right) further 
degrades the perturbative result. 

expansion seems to hold as an asymptotic series. The inclusion of higher order corrections 
degrades the result and the closest we get from the exact numerical result is by keeping 
0(2) corrections (see figure 3). However, even though the perturbative result fails for 
large a, it is consistently larger than the numerical one, which is supposed to be the true 
minimum. The lowest order correction shows that for small a, the distortion is suppressed 
for large number of extra dimensions. Furthermore, for incidence along the vertical axis, 
the distortion is independent of £ (this follows trivially from the azimuthal symmetry). 
For incidence along the equatorial plane (or intermediate angles), the disc of absorption 
is distorted into an oval, indicating the orientation of the angular momentum of the black 
hole. The exact numerical results 2 in figure 4 confirm the non-perturbative validity of 
these qualitative features. The only exception is the conclusion regarding the suppression 
of the distortion with n (actually the opposite occurs for a* > 1, see right plot in 4 where 

2 These were obtained by running a code in Maple which was written using the Maple inbuilt minimisation 
routine. 
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Figure 4. The plots show the variation of the absorptive disc by varying one parameter with all 
the others fixed. 

the distortion is greater for larger n when a = 4). This figure also shows that for larger 
rotation parameter, the distortion becomes much stronger. 

To better understand the contributions to the fluxes at high energy, it is instructive to 
look into an approximation based on the geometrical cross sections obtained from the discs 
above. In [59], Unruh proved that the absorption cross section for a plane wave, incident 
on a Schwarzschild black hole from infinity is 3 



This suggests the interpretation of Oj^ m as the contribution from a partial absorption cross 
section for a wave with angular momentum quantum numbers {j, m} . These cross sections 
are directly related to the transmission factors so they are usually named greybody factors, 
because they are responsible for distorting the black body spectrum of Hawking radiation. 
In this non-rotating limit the flux becomes 



so the Planckian term factors out of the sum. In the limit oo — > +oo we would expect 
an incident plane wave to be well described by a beam of classical particles. Then the 
absorption cross section is simply the area of the absorptive disc at infinity (which in the 
non-rotating case is a circular disc) . This type of approximation was noted by DeWitt [60] , 
who replaced the transmission coefficient by a theta function cutting the j-sum in (2.18) 
at the maximum angular momenta allowed by the absorptive disc radius (j = b max uj). The 
success of his approximation, reinforces the interpretation of a as the total cross section 
for a classical beam of particles. 

Expression (2.18) is more complicated for a rotating black hole (equation (2.5)) because 
the Planckian factor depends on m and can not be taken out of the sum. However, it is 
still possible to prove that for a wave incident at an angle 9, the cross section is (we are 

3 We are using horizon radius units. 
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Averaged geometrical cross-section 




a 



Figure 5. The plot shows a perfect agreement between the averaged geometrical cross section 
(solid lines) and the asymptotic value of the sum of the greybody factors (the points). We show the 
cases from n = 1 (red upper line) to n = 6 (bottom line). The curves are organise in order from the 
top to the bottom. We also show the n = 1 case which could be relevant for the Randall-Sundrum 
model [5, 6], though we do not consider it in the remainder. 



using the scalar case) [61] 



< e ) = ^T Yl \Sk(c, cos d)\ 2 T Cik = °<u,mi6) (2-19) 

k={j,m} k={j,m} 

if we average over the solid angle we get that 

4^)= E ~^ k= E ( 2 - 2 °) 

k={j,m} k={j,m} 

So for the rotating case, the same relation between partial wave cross-sections an trans- 
mission factors exists. But in the high energy limit, assuming the geometric description 
is good, we know how to compute the absorption cross-sections for a given angle of inci- 
dence by using the discs obtained above. Furthermore, note that at high energies we would 
expect this result to be independent of the spin of the particle. In figure 5 we compare 
the geometrical result obtained by numerically integrating the absorptive discs, with the 
asymptotic value of the sum over transmission factors computed for the scalar field. We 
find an excellent agreement between the points taken from Table I of [55] (where a wave 
scattering calculation was done) and our geometrical calculation. 

Furthermore, one can extend DeWitt's argument to compute the Hawking flux. If we 
go back to the scalar angular flux before integration over the solid angle d£l we have 

dN 1 ^ 1 



dtdcudtl 4tt 2 ^-^ exp^Cj/Tn) 
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Since we know the partial absorption cross sections as a function of the angle by using the 
geometrical discs, we can find a high energy approximation for the spectrum by cutting off 
the sums according to the allowed regions on the (b, Q plane. A particularly interesting 
limit is when we approximate the exponential in equation (2.5) as 
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which holds for a small or large or for small enough j. Furthermore j ma x is large, when 
oj — > +oo, so for most of the modes contributing to the the sum in (2.5) this approximation 
should work. This implies a similar factorisation in the high energy limit for rotating black 
holes: 



dN 



OJ 



dtdoj 2vr 2 exp(oj/T H ) - 1° ' 

The approximations obtained in this section provide some interesting information 
about the Hawking spectrum in the high energy limit (where it should be independent 
of the spin of the particle). As expected, the discs for incidence along the angular momen- 
tum axis are circular, in agreement with the symmetry of the background. The angular 
asymmetry as we vary 8 becomes manifest in the shape of the discs. As verified in (2.20) 
and (2.24) there is a direct correspondence between the absorptive discs and the Hawking 
spectrum at high energies. Another interesting property that we can see in the middle 
plot of figure 4 is an enhancement of the Hawking emission in the high energy limit, as we 
increase the rotation parameter of the background (another known generic property of the 
Hawking spectrum), since the discs grow in size. Some of these properties will appear in 
the next section. 



2.2 Exact Hawking fluxes 

The angular spheroidal functions S S (aoj, cos 8) are the main ingredient inducing non- 
uniform angular distributions which we shall investigate. They are obtained by solving 
the spheroidal wave equation as detailed in previous studies [22, 24-27, 30-32, 56]. To 
quickly recall the relevant properties for our analysis, we will simply use the data avail- 
able in [62] for the transmission factors and the method in [36] to evaluate the spheroidal 
functions and obtain the angular fluxes for different brane degrees of freedom. 

Some important properties of the angular distributions follow from the observation of 
the differential number fluxes in figure 6. It is well known that in general rotation tends 
to make the spheroidal functions more axial whereas the combination of the transmission 
factor and Planckian factor favour m = j modes at large rotation to spin down the black 
hole. The competition between the angular functions and the latter results in more axial 
angular distributions at low energies and more equatorial distributions at high energies 
when rotation is on. This is seen from the energy dependence of the angular profiles shown 
in figure 6. 



- 10 - 




Figure 6. Angular dependence of the flux for states of positive helicity h = — s = 0, 1/2, 1, non-zero 
rotation a = 1 and two different n = 2,6. The negative helicity plots arc obtained by the reflection 
cos 9 — > — cos 9 

These effects are particularly large for vector bosons, which are more likely to be 
emitted close to the rotation axis at low energies, whereas high energy vector bosons are 
more likely to be emitted in the equatorial plane (the same happens for fermions). In 
figure 6 we see that particles with a single helicity will be emitted asymmetrically by 
a rotating black hole [63]. For example, if the black hole angular momentum vector is 
pointing north, then positive(negative) helicity states will be preferentially emitted in the 
northern (southern) hemisphere. As suggested before in [36], this might lead to angular 
asymmetries in the decays of unstable W and Z vector bosons which can be in two possible 
transverse polarisations states. 

3 A parton level study of angular correlations 

The remainder of this article is dedicated to studying the feasibility of observing these 
angular asymmetries in a more realistic model for the decay of TeV gravity black holes at 
colliders. We use the recently upgraded CHARYBDIS2 event generator which includes in full 
the effects of rotation for brane degrees of freedom, including the helicity effects described 
in the previous section. We have produced some samples to illustrate the angular effects in 
question. These effects should be fairly independent of the details of the sample, given that 
rotation is expected to be present in general, and the qualitative features of the angular 
fluxes described in the previous section are fairly independent of n. 

3.1 Angular momentum reconstruction 

The problem of determining the angular momentum axis and /or magnitude is much more 
difficult than reconstructing the mass of a rotating black hole. This is because the angular 
momentum axis, which controls the angular distributions, evolves during the evaporation, 
changing direction and magnitude. Furthermore, all the decay products in the laboratory 
frame are in general boosted with respect to the black hole centre of mass frame, and the 
black hole recoils between each emission. Nevertheless, we may hope to see some of the 
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effects by boosting each event back to the partonic centre of mass frame. In principle we 
can determine the latter reliably if there is little missing energy emitted in the evaporation. 
The only missing energy in the SM comes from neutrinos and they represent a small part 
of the total number of degrees of freedom. However missing energy from neutrinos can 
arise from the secondary decay of SM heavy particles. Nevertheless the fraction of missing 
energy from the evaporation should not be so large (for most of the events) as to degrade 
the results much more compared to other factors such as the recoil. Note however that 
including gravitons in the evaporation may degrade the reconstruction further, though in 
principle, by measuring all the proton remnant pieces we may be able to obtain a good 
measure of the partonic centre of mass frame. Alternatively to reduce this effect we may cut 
on events with a small amount of missing energy as it has been done in previous analyses 
for the mass reconstruction (see [36] and references therein). 

We start our study by considering some distributions which are not accessible observ- 
ables at a collider experiment, but help to understand the angular correlations involved. 
The aim is to understand how the momenta of the particles emitted correlate with the true 
angular momentum axis by using knowledge of the black hole decay history. 

The first distribution we study is the angle of a particle with a given spin with the 
true initial angular momentum in the centre of mass frame of the initial black hole 



The first plot in figure 7 shows a spin dependent behaviour as expected. The sample 
used to produced the plots used all the defaults values except for M JL0ST= . FALSE . . This 
produces black hole events with larger angular momentum and mass. This is not essential 
since we could use a cut on the visible invariant mass for the event to select heavier black 
holes if M JL0ST= . TRUE . , it is simple advantageous to generate samples faster. We know 
from the theoretical plots in figure 6 that scalars and fermions tend to be more equatorial 
(though fermions at low energies also have a small axial peak) and vector bosons are very 
axial at low energies. This is consistent with the larger probability for vector bosons at 
larger |cos6>| in the top left plot (though the effect is not very large). We can improve 
this correlation by selecting particles in particular ranges of energy for each initial black 
hole. For example if our cut requires high energy particles for all spins we would expect an 
equatorial correlation, whereas at low energies we expect a flatter distribution for scalars 
and fermions and an axial distribution for vector particles. This is confirmed in the two 
top right plots and the bottom plots of figure 7 where we have chosen ranges of energy in 
units of the horizon radius of the initial black hole. This result suggests using soft vector 
particles as a guess for the axis, to plot angular distributions. This will fail a considerable 
part of the time and smear out the true correlation. 

High energy particles tend to be emitted perpendicularly to the initial black hole 
angular momentum. This is seen in the bottom right plot of figure 7 where the correlation 
is stronger. If in addition we assume that the direction of the angular momentum vector 
is perpendicular to the direction of the black hole momentum (which is true in the limit 
where the angular momentum does not recoil during the production) we obtain another 
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Figure 7. Distribution of angles of particles with various spins with the angular momentum axis. 
The top left plot shows the distributions for cos 9 for various spins for an n — 3 sample with black 
hole masses above 5 TeV. The other plots show the distribution for various ranges of energies as 
indicated in the title of each plot. Note that the error bars for s = are due to the small number 
of pure scalar degrees of freedom in the SM (only the Higgs particle) . 
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Figure 8. Axis reconstruction: Using the hardest emission (left) and the eigenvector of the smallest 
eigenvalue of the sphericity tensor (right). 



guess for the axis. The left plot in figure 8 shows that this method works better. However 
it relies on the assumption that the initial angular momentum is perpendicular to the black 
hole momentum. 

A third method to estimate the angular momentum axis is to consider the shape of 
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the event, i.e. to use all the momenta in the decay. For a rotating black hole we expect 
most of the particles to be emitted equatorially (since only low energy vector bosons are 
axial). Thus the event should have a disc like distribution of momenta indicating the 
orientation of the axis. This axis should minimise the amount of momentum projected 
along its direction. Another advantage of this reasoning is that it gives lower weight to 
low energy particles which we want to eliminate since low energy vector bosons are more 
axial, and low energy scalars and fermions are more uniform. If we denote the direction of 
the angular momentum by n, and use projections of momenta squared, then we want to 
minimise 

min n £i(Pi • n) 2 _ min n £\ J2 af s n °P? P? n/? 



• \^ a ( Ei PfP^ 
= mm 2 J n 



EJp 




minrn^n" (3.2) 

n ^ — J 



n 

a/3 

where we are use Greek letters for spatial indices and the index i runs over all particles in the 
event. We have defined the sphericity tensor S a @ as usual [64]. The sphericity tensor has the 
properties that all eigenvalues are non-negative and their sum is one. In the eigenbasis it is 
clear that the direction which minimises the quantity in (3.2) is the eigenvector associated 
with the smallest eigenvalue. In the right plot of figure 8 the distribution for the angle 
between the guessed axis (the eigenvector associated with the smallest eigenvalue) and the 
angular momentum axis is shown. This method is not as good as the the left plot but it 
has the advantage of not relying on the assumption that the angular momentum is on the 
plane transverse to the collision axis. 

3.2 Angular correlators 

An alternative to reconstructing the black hole angular momentum to study angular dis- 
tributions, is to explore the property that polarised angular distributions are strongly 
dependent on the helicity of the particle. For example from figure 6 we know that there is 
a strong preference for vector bosons to be emitted in different hemispheres. This motivates 
defining angular correlators of the form (in the frame of the initial black hole) 

= 1 Ti 7 (^'^) 
IPillPjl 

which are cosines of angles between particles % and j. So for the case of particles with the 
same helicity, we would expect the distribution to be higher at Xij ~ 1 and reduced at 
Xij ~ — 1, and the opposite to happen for particles with opposite helicities. Figure 9 shows 
the expected behaviour for pairs of emissions in a fixed black hole background (no recoil). 
The probability density function used to determine the distribution of (3.3) is derived in 
appendix B. The effect at fixed black hole parameters grows quickly with a*, especially for 
particles of helicity h = ±1. The bottom right plot shows that there is some variation with 
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p ) fol " = 4 p (* 1 > 1 ) for n = 4 p ^°>°) for n = 4 




Figure 9. Angular correlators for an eternal black hole. The top plots show the probability density 
functions for the angular correlator between particles of same helicity for fermions, vector bosons 
and scalars (left to right), for a range of a* values. The two left bottom plots are similar but 
between particles with opposite helicities. The bottom right plot shows the dependence with n for 
a* = 1 and two vector bosons with the same helicity. 



n though not very strong (the curves are qualitatively similar), as claimed in the beginning 
of this section. 

For an evolving black hole, we expect these effects to get smeared, due to the mo- 
mentum and angular momentum recoil. Again the best procedure is to compute similar 
quantities in the rest frame of the initial black hole assuming a small amount of missing 
energy during the evaporation. 

Figure 10 shows the correlators defined in (3.3) for various helicity combinations. The 
top left plot shows the distribution for same helicity (red) and opposite helicity (blue) 
fermions (solid lines) and vector particles (dashed lines) and the non-rotating case (black) 
for comparison. The black curve is not constant due to the recoil of the black hole be- 
tween emissions. The recoil tends to make different pairs of emissions more back to back 
(especially subsequent emissions) which causes the rise at Xij ~ — 1 and the fall off at 
XiJ ~ 1. 

We can see (particularly for vector particles) that the asymmetry predicted in figure 9 
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Figure 10. Parton level angular correlators evaluated in the frame of the initial black hole. The 
sample used was the same as in figure 7. The top left plot shows correlators between any two final 
state fermions or any two final state vector bosons for some helicity combinations. The correlator 
for any two particles when rotation is off is shown for comparison. The top centre and top right 
plots show the same distributions using various pairs of particles according to their ordering during 
the evaporation. The central row of plots shows the same distributions using specific intervals of 
energy (compared to the horizon radius for the initial black hole - rno)- Twice the bin size was 
used to compensate for the lower statistics. The bottom row of plots contains the same plots as 
the top row, but with a different recoil model REC0IL=2. 



persists, with an enhancement at Xjj = — 1 for opposite helicity correlators, compared to 
the same helicity correlators. It is also clear that the effect of the recoil between emissions 
is larger in the rotating case. This is because the spectrum is harder, hence the relative 
boost when the black hole recoils is also larger. Additionally, in the rotating case, we have 
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a forward peak at Xij = 1. This is due to the possibility of having an accumulated boost 
in a certain direction for the particles emitted later in the decay, so pairs of those particles 
will tend to be more collinear. Since in the rotating case the multiplicity of the event 
is reduced and the spectrum is harder, this effect tends to be important (practically all 
events end up with a large accumulated boost). This is verified in the top centre and top 
right plots where for the first pair in the beginning of the evaporation (top centre plot) no 
forward peak is observed, whereas for later consecutive pairs, the forward peak tends to 
increase rapidly. For pairs which are increasingly separated in the order of emission, the 
forward peak also disappears (top right plot). The reason why no forward peak appears 
in the non-rotating case can be justified as follows. For events with multiplicity N, the 
number of consecutive pairs which can contribute to a forward peak is iV — 2, whereas the 
total number of pairs is N(N — 1) /2. So the fraction of pairs contributing to a forward peak 
is at most 2{N — 2)/(N(N — 1)). Typically, when rotation is turned off, the multiplicity 
increases from 5~8tol0~15so the fraction of pairs contributing to the peak is reduced 
roughly by a factor of ~ 2. Furthermore, a larger multiplicity means a smaller magnitude 
for the boost in the recoil by another factor of ~ 3 (since the same energy is distributed 
among more particles which are softer). So overall, we have a suppression factor of at least 
~ 6. This explains the absence of the forward peak for the non-rotating sample. 

The central row of plots shows that by selecting particular ranges of energy we recover 
the strong asymmetry for low energy vector bosons (left and centre plots). For high energy 
particles (right), all helicities are equivalent. This is because at high energies all the angular 
spectra become equatorial regardless of the spin. 

The bottom row plots are the same as in the first row, except that REC0IL=2 was used 
(see section 3.3 of [36]). Similar conclusions are obtained with this option (particularly for 
the ranges of energy selected in the middle row of plots). The only difference is that the 
forward peak due to the recoil is somewhat sharper and the backward effect is smaller (see 
explanation in appendix C). 

4 Conclusions 

In this article, we have performed the first parton level study of angular asymmetries in 
semi-classical black hole events at the LHC. 

In the first part, we have done the first complete analysis of the geometrical optics discs 
for our background, which provides information on the Hawking fluxes at high energies, 
and a picture of the background as seen by an asymptotic observer. The calculation was 
shown to match with some properties of the full Hawking fluxes so it provides an alternative 
to model this asymptotic regime. 

Phenomenologically, we have devised several strategies to extract angular information 
from the events. In the first type (axis reconstruction) though the reconstructed axis 
correlates with the true axis, there is a large fraction of wrong reconstructions due to recoil. 
This means that an angular distribution of the momentum of a particle with respect to such 
axis, is likely not to present a strong correlation. After observing this, we have turned to 
angular correlators between pairs of particles, which do not depend on axis reconstruction 
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and explore much better the angular asymmetries of the events. We have devised a set of 
energy cuts which are in principle easily generalised to a realistic experimental situation 
by using the energy scale of the process (which determines rno)- 

We should note that the plots in figure 10 were produced assuming knowledge of 
the helicities of the outgoing particles. This is usually not a direct observable at hadron 
colliders. However the decay modes of the the W and Z bosons are dependent on the helicity 
of the intermediate states. For example in the decay W~ —> i~vi the charged lepton tends 
to be collinear with the W~ for negative helicity, and anti-collinear for positive helicity. 
For Z decays a similar argument holds. The only disadvantage of W decays is the neutrino 
which makes re-construction harder, whereas for Z decays the disadvantage is the mixture 
of left-handed and right-handed couplings. Finally, gluons, being vector particles should 
also exhibit similar asymmetries, so cutting on lower energy jets could potentially produce 
similar asymmetries, even without helicity knowledge. 

To conclude, our analysis shows that there is an effect which in principle can be ex- 
tracted from the data at the experimental level should black hole candidates be found in 
the future runs at the LHC. A full detector level analysis, which is beyond the scope of 
this work, will tell us how to extract in detail these interesting effects. 
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Appendices 

A Critical impact parameters 

In this appendix we describe how to obtain the critical impact parameter discussed in the 
main text. When Q = 

A = r n+2 + r n a 2 (l - cos 2 C sin 2 0) - (1 + a 2 )r 
B = — a cos C sini? (l + a 2 ) r 

C = r n+4 (l - v 2 ) + r n+2 a 2 (l - v 2 + cos 2 ■&) + r n a 4 cos 2 + 

+rV(l + a 2 ) + ra 2 (l + a 2 )sin 2 i? (A.l) 

where C > since v < 1. For a fixed b, it can be shown that A has exactly one zero in the 
domain r > so it changes sign only once 4 . So A starts from positive values at infinity, 
decreases, goes through zero and takes negative values inwards. Since the signs of B and 
C are fixed we only need to analyse two regions. 

4 This is done by looking at the sign of the function at infinitesimal r and at infinity, using continuity 
and the positivity of d 2 A/dr 2 > for r > 0. 
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Let's start with the outermost region. We want to determine the impact parameters 
for which the particle is not able to penetrate completely through the region of A > 0, thus 
being scattered back at some point, or equivalently the range of parameters b m i n < b < 

u max 

for which absorption is guaranteed. Then it is certain that any particle with impact 
parameters within this range will reach the second innermost region. For any solution 
of the radial equation we must ensure that 1Z > to allow the particle through the region. 
As a function of b, 1Z has two zeros 

b±= B ±V^TW (A2) 

Since, AC is positive, there is a negative and a positive root. Thus 1Z is a parabola with 
a maximum. 1Z is positive only for < b < 6_|_. If we take b > b max = min r {6+}, where 
the minimisation is over the radial region we are considering, there will always be a point 
where 1Z goes through zero and changes sign. This means that the particle is scattered 
back at that point. Similarly, in the complementary case, there will never be such a zero 
so the particle reaches the second interior radial region. Thus we have obtained the first 
upper bound. 

Regarding the remaining interior radial region the situation is a bit more complicated. 
In that case, A < so the roots are 



i± ^ Z B±v*Em . (A . 3) 

For B > (or equivalently b z < 0, see (2.11)), the previous equation has no positive real 
roots, so 1Z > and any particle reaching the region is absorbed. Thus, for this sign of B, 
the absorptive disc is defined by the interval obtained in the outer region. For B < 0, we 
can have two real positive roots. Those are 



6t = im ^Em (A „ 

so 7Z is a parabola with a minimum and two zeros. If b takes any value below min r {6_} 
(again the minimisation is in the radial region we are considering), there is no zero, because 
any r has associated critical 6's which are necessarily larger. The same occurs for b above 
max r {& + }, since all the Vs that allow a zero are smaller. Conversely, for b between the 
last two values, there is always a point where 1Z goes through zero and becomes negative, 
so the particle is scattered back. Therefore, the impact parameters for absorption must 
be smaller than min r {6-} or larger than max r {6+}. However, note that at r = ta-, (r^ 
defined such that A{ta) = 0) b + diverges, so there is no finite b > max r {&+} — > +oo. 

A.l Recurrence relations 

Since r c , b c are the minimiser and minimum when a / 0, r c must obey 

| =0, (A.5) 
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so that 



b \r=r c 



(A.6) 



Note that from now on, b is given by the expression to be minimised in (2.13) and all 
other parameters f , n are omitted. Expanding (A. 5) in powers of a and imposing each 
coefficient to vanish independently, we obtain conditions on the r, p coefficients (see equa- 
tion (2.16)) 



d m 


( db 






da m 


. dr 
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a=0 



. 



(A.7) 



One can show that the total derivatives can be expanded in the following manner 
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By successively iterating (A. 8), one can find a general expression which is proved by in- 
duction 



d m 
da™ 



m—l 

E 

g=0 



(m — 1)! / dr, 



- 



(m — 1 — q)\q\ \da J da m q dr. 



r + 



Y^ 1 d m - q r c 



dr c 



Qq+1 



3 dam ~ q t^o ( q - k)[kl ^ da ' da q - k dr k + l 
Using (A. 9) in (A.7), relabelling r c — > r in the partial derivatives, and defining 

d i+j b 



(A.9) 



da l dri 



(A.10) 



a=0,r=ro 

gives the recursion relations (m > 1) for all the corrections 

m— 1 
9=1 

The first correction is 



m—l 



r\p m - q , q+1 + r m _ q ^ r\ ( Q \ f3, 
k=o ^ ' 

2 cos C sin "& 



q-k,k+2 



V(n+l)(n + 3) 

Furthermore, a similar expansion for b c can be found by using (A.9). Define 

oo , 

ml 



m=0 



(A.ll) 



(A.12) 



(A.13) 



where 



d m+l b c 




' d m / 


da m+1 


a=0 


da m \ 



' d dr c d 
da da dr c 



J a=0 



d m db c 
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(A.14) 
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The last step follows from (A. 7). The result is (m > 1) 



h = /3i,o 

m— 1 
9=0 



m — 1 



T\@m-\-l—q,q Tm—q 



ty+l-fc,fc+l 



(A.15) 



The first correction is actually independent of the perturbation in r 

2 cos C sin $ 



6i 



n + 1 



(A.16) 



B Eternal black hole angular correlators 



In this section we obtain the angular correlators on an eternal black hole background, used 
in the main text. We define the probability density function for the correlator Xij (up to 
a normalisation constant) 

p(xij) oc / dxidxjd(t)id(t)jPi{xi)pi(xj)5 (x itj - yjl - x\ Jl - x 2 - cos((^ - <j>j) - XiX. 



where we have defined the spatial momenta of particle i (or j) 



1 - xfcosfa, a/1 - xfsmfii^i |p 



and 



+oo „ 

pi( x i) = J2 



dbJ 

K=0 



l S k (au},Xi)\ z 



exp(bjr H /T H ) ± 1 



(B.l) 



(B.2) 



(B.3) 



is the probability density of having a particle of type i (with any energy) emitted with 
direction Xi with respect to the angular momentum axis. Due to the azimuthal symmetry, 
the distribution is uniform in the 4>i direction. This can be written in a more convenient 
form by using the definition 



uir H 





1 



to obtain 



+oo „ +00 

Pi(xi) = yz 

K=0 J ° 



exp(x/Tn) ± 1 ~ 
dC h>K ^ D (u)r H ) 



Tr\x, a*) 



doj 



\S k (auj,Xi)\ . 



(B.4) 



(B.5) 



Then making the change of variable y = fx{u) = Ch,K,a,,D{'^i'H)/Ch,K,a,,D{+ ( X') 

-oo r i 
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where /( 1 ) is the inverse function (not 1//). Now if we define 
Ci{xi) = / dxpi(x) 



i 

= *E C ^„f)(+») dxISfcCa/^^)^)! 2 (B.7) 
Jo K=Q J-i 

=> = • (B.8) 

Going back to equation (B.l) we perform the changes of variables 

$ = & + <f>j (B.9) 

= 4>i - <f>j (B.10) 

= 9i(xi) = Ci(xi)/Ci(l) (B.ll) 

to obtain 

-1 fl fin f -2n 



p(xij) oc J J J J dwi dwj d<& dcj) 5 ^Xij — \j \ — x 2 — x 2 cos 4> — XiXj^j 

p(Xij) = j j j dWi ^ & (^ X i,j - \Jl - xfijl - xj COS (f> - XiC 



1 fl f2lT 

(B.12) 

where x$ = g^^iwi) and we have normalised the distribution. The histogram for p(xij) 
is obtained by generating the phase space Wi, Wj, <fi uniformly and adding a unit weight to 
the bin for the corresponding 



1 — xfyl — x 2 cos 4> + XiXj . (B.13) 
C The effect of recoil 

In this section we explain the backward and forward peaks found in figure 10, using the 
kinematics of the two recoiling options in CHARYBDIS2. Using the description in [36], it can 
be shown that the mass reduction for a particle with a small energy E <C M (i.e. in the 
beginning of the evaporation) in the frame of the initial black hole is 



M 



final 



M - E ( 1 + -4t 1 + • • • , RECOIL = 1 
M — E , RECOIL = 2 



This explains the larger backward effect (at x = —1) in the top plots when REC0IL=1, 
since in that case the black hole mass reduction is a bit larger in the beginning of the 
evaporation, enhancing the effect of the back to back recoil. The sharper forward peak 
when REC0IL=2 can be explained by expressing E in terms of the selected Hawking energy 
u (neglecting particle mass) 

U3 UJ 



E 
M 



r £ [0, i] , RECOIL = 1 
M M 1 21 

/ 1 \ • (C.2) 
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BH frame energy vs Hawking energy 
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Figure 11. Energy of the particle emitted: The energy of the particle (E/M) in the frame of the 
black hole is shown as a function of the energy selected from the Hawking spectrum (lu/M). 

It is easy to see (figure 11) that at high energies E close to the kinematic limit (which are 
more important in the last part of the evaporation), the range of energies u contributing 
to a range of energies E is always much wider for REC0IL=2. So there will be more 
hard particles selected at the end of the evaporation for the latter, contributing to the 
accumulated boost and hence the forward peak. 
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